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t^. ■ Abstract 

Macroscopic derivation of the entrainment in a supersolid cylinder induced by a surface 
elastic wave (of small amplitude) propagating along the flexible interface is conducted by 
Qh considering the nonlinear coupling between the interface and the rarefaction effect. We ob- 

tain the critical bounds for zero- volume-flow-rate states corresponding to specific rarefaction 
measure and wave number which is relevant to the rather small critical velocity or disap- 
pearance of supersolid flows reported by Rittner and Reppy. 
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In 1969, it was conjectured by Andreev and Lifshitz 1 that at zero temperature, delocalized de- 
fects may exist in a quantum solid, as a result of which the number of sites of an ideal crystal 
lattice may not coincide with the total number of particles. Originally, this conjecture was pro- 
posed for three dimensional quantum solids made of atoms ( 3 He, 4 He, • • •) which do not interact 
via Coulomb repulsion. The proposed supersolid phase is believed to occur due to the quan- 
tum behavior of point defects, namely vacancies and interstitials, in this crystal of bosons 2-3 . 
Researchers have found that a small lattice model does not exhibit the mesoscopic signature of 
£>V an intermediate phase separating the solid from the liquid, where the solid and the fluid would 

coexist 4 . Such a vacancy-solid phase was indeed suggested 1 by Andreev and Lifshitz if the zero 
point motions of certain defects become sufficient to form waves propagating inside the solid. 
Castaing and Nozieres have later considered 5 such a possibility for spin polarized 3 He. The 
statistics of the defects depend on their nature. For simple vacancies in the crystal, their statis- 
tics is given by the statistics of the particles out of which the solid is made. If the defects are 
bosons, they may form a condensate, giving rise to a superfluid coexisting with the solid. This 
supersolid phase is discussed in certain bosonic models 6 . If the defects are fermions, they may 
form a Fermi liquid 7 coexisting with the solid, such that the system is neither a solid, nor a 
liquid. Two kinds of motion are possible in it; one possesses the properties of motion in an 
elastic solid, the second possesses the properties of motion in a liquid. This interesting issue 
motivates our present study. 

Early theoretical work by Andreev and Lifshitz 1 and Chester 2 showed that solids may feature 
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a Bose-Einstein condensate of vacancies (or interstitial atoms) and thus possess superfluid (SF) 
properties. Quite recently one description of the quantum solid is as a density wave that has 
formed in the quantum fluid 8-9 . The periodicity of this density wave need not match precisely 
to the particle density, so that the ground state may be incommensurate, with unequal densi- 
ties of vacancies and interstitials. Whether or not the same is true for quantum fluctuations 
is not clear at this point. We noticed that previous theories imply a corresponding vacancy 
contribution to the specific heat that is as large as the phonon contribution near 1 Kelvin 3,10 . 
Based on these considerations or phenomenological approaches, assuming the existance of small- 
amplitude density waves along the deformable boundaries, in this letter, we shall demonstrate 
that wavy flexible interface (between atoms and free vacancies or defects) or highly-pressured 
environments 8 can produce elastically deformed interface or peristaltic motion will induce time- 
averaged transport in a Andreev-Lifshitz supersolid 1 ' 2 . 

Theoretical studies of interphase nonlocal transport phenomena which appear as a result of a 
different type of nonequilibrium representing propagation of a surface elastic wave have been 
performed before 11-12 . These are relevant to particles flowing along deformable elastic slabs with 
the dominated parameter being the Knudsen number (Kn = mean-free-path/L^, mean-free-path 
(mfp) is the mean free path (of the particles) which is thus temperature dependent 13 ' 14 , is pro- 
portional to the characteristic distance between two boundaryies) 13-15 . The role of the Knudsen 
number is similar to that of the Navier slip 15 parameter N s (= /jS/Ld] S is a proportionality 
constant as u s = St, t : the shear stress of the bulk velocity; u s : the dimensional slip velocity; 
for a no-slip case, S = 0, but for a no-stress condition. S = oo, \i is the viscosity). 
We shall choose a periodic domain to simplify our mathematical treatments. The deformable 
(elastic) interface is presumed. We adopt the macroscopic or hydro dynamical approach and 
simplify the original system of equations (related to the momentum and mass transport) to 
one single higher-order quasi-linear partial differential equation in terms of the unknown stream 
function. In this study, as the temperature is rather low and the phase is related to the super- 
solid (there might be weakly friction or shearing dissipation in-between) we shall assume that 
the governing equations are the incompressible Navier-Stokes equations which will be associated 
with the microscopically slip velocity boundary conditions along the interfaces 13-15 (cf. Ref. 
13 for the quantum slip case). To consider the originally quiescent environment for simplicity, 
due to the difficulty in solving a fourth-order quasi-linear complex ordinary differential equation 
(when the wavy boundary condition are imposed), we can finally get an analytically perturbed 
solution and calculate those physical quantities, like, time-averaged transport or entrainment, 
critical forcing corresponding to the freezed or zero-volume-flow-rate states. The latter might 
be relevant to those reported in Refs. [16-18] for very low flow rates 16 or disappearance of su- 
persolidity through an annealing of the solid helium sample 17 . Note that, as also reported in 
Ref. 18, their results suggested that grain boundaries (GBs) are superfluid, so that 4 He crystals 
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of medium quality are supersolid at the liquid/solid equilibrium pressure (P m ); that is, mass 
transport through them without dissipation is possible. 

2 Formulations 

To escape from the difficulties in treating many-body or many-particle problems together with 
the scattering with a elastic boundary, we adopt the hydrodynamical approach but use the 
microscopic quantum slip boundary condition which takes into account the mean free path 13 of 
the dilute molecular gas corresponding to the nonzero slip velocity 13 along the interface-wall. 
We consider a circular cylindrical tube (that of the vortex core) of uniform radius filled with 
a homogeneous rarefied gas (Newtonian viscous fluid). The wall of the tube (or the interface 
between the inner vortex core and the outside part of it) is not absolutely rigid, on which 
is imposed axisymmetric travelling sinusoidal waves of moderate amplitude a (z is the axial 
coordinate in the wave propagation direction). The radial displacement from the mean position 
of the wall or the interface (r = r w ) is thus presumed to be r/, where r] = acos[27r(z — ct)/X\, 
A is the wave length, and c the wave speed (cf. Fig. 1). Axisymmetric motion is assumed 
with r measured in the direction normal to the mean position of the wall, u, v are the velocity 
components in the z- and r-directions, respectively. 



a = e r„ 



T 



rf(z,t) = a cos x ( z — c *) 



wave speed c 



Fig. 1 Schematic diagram of the wavy motion of the interface or tube wall. e<l 

Firstly, it is necessary to simplify these equations by introducing nondimensional variables. We 
have a characteristic velocity c and three characteristic lengths a, A, and r w . Thus, the following 
variables based on c and r w could be introduced : 

, r , z ,'/,'',//, tp , ct , p 
r = — , z = — , u =-, v = -, rj = — , V = — n"> * = — i P = — 2> 

w ' w ° ° ' w c ' w ' w 

where is the stream function, p,p are the density and pressure of the fluid. The amplitude 
ratio e (presumed to be rather small), the wave number a, and the Reynolds number Re are 
defined by 

a 2nr w cr w 

e = — , a = — - — , He = . 

r w A v 
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We shall seek a solution in the form of a series in the parameter e : 

,1,1,2, dp dp dp 2 dp 

ip = ip + eV'i + e V2 + • • • , 7^ = (7^)0 + e (^)i + e (^)2 + • • • . 

with v = (dtp/dz)/r, u = —(dip/dr)/r. 

The r- and z-momentum equations and the equation of continuity 19,20 could be in terms of the 
stream function ip if the pressure (p) term is eliminated. The final governing equation is 

» *V + 'M^a - ?vv + % - = ^v>, = |i + |i - if (1) 

ar r r r z r fie dz z or z r or 

and subscripts indicate the partial differentiation. Thus, we have 

—V Vo + V Vo r - - V Vo + —5- V Vo z = -75- V ip , 2 

at r r r z r Re 



dt 



ry /y /y A ly 

— v Vi r - -v ^1 + — v v>i z - D 



[V Vi r - -V ipi + -5-] V ip lz = — V Vi, (3) 



_v ^ + —[V v 2r - -v V2 + ^2-] - —v v 2z + —[v v>o r - -v Vo + 



" — V V'o. + — [V - - V Vi + ^2-] - —V Vl, = ^2, (4) 



and other higher order forms. The fluid is subjected to boundary conditions imposed by the 
symmetric motion of the walls and the non-zero slip velocity 13 ' 14 : u = -Kn du/dr, v = drj/dt 
at r = (1 + 77). The boundary conditions may be expanded in powers of r\ and then e : 

— Id 1 d d 2 

[( ~Tfr^° + ^ + ^ + ' ' ■)]lr=l+ e ooea(z-t) = " Kn [(^^ " ^)(V>0 + 

6^1 + e 2 ^ 2 + • • -)]lr=l+e cos a(z-t), (5) 

ip0z\i + e[cosa(z - i)^ozr|i + ipu\i] + e 2 [ 3^X11 cos 2 a ( z - t) + cosa(z - i)V'Ur|i+ 

^2z|i] + ••• = £« sin a(z — t) + e 2 a cos a(z — t) sina(z — t). (6) 

Equations above, together with the condition of symmetry and a uniform pressure-gradient in 
the z-direction, ((9p/9z)o=constant, yield : 

V>o = ^o[r 2 -y], K = f(f z ) , (7) 

Y>! = 4>{r)e ia ^ + 4>*{r)e~ ia ^ + fo(r), (8) 



where the asterisk denotes the complex conjugate. A substitution of ipi into equation (3) yields 

{-T2~ a2 -^r + iaRe[l + 2Ko(l-r 2 )}}(-^-a 2 - -i-)<f> = 0, V 2 V 2 </> = 0. (9) 
dr z rdr dr z rdr 
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The boundary conditions are 

<Mi = 2K , Or |i = 4K o Kn, <t>(l) = ~. (10) 

The equations for </>* are conjugate to the foregoing and need not be written down. 
Similarly, with 

^ = D{r) + E(r)e i2aiz - t) + E*(r) e - i2a ^ + G(r)e ia ^ + G* (r)e~ ia(z ~ t \ (11) 



we have 



d 2 d d 2 d d 1 1 ?'o-/?p 

rr rrrr , \{4>€ ~ W)] - iaRe{^[ ^f^ ] ~ @lZJLlM } . (12) 

(TT " 4 « 2 " A~? E = -&<*Re[l + 2K (1 - r 2 )](^ - 4a 2 - -^-)E+ 
dr z rdr dr A rdr 

iaRe r , , ,, 3,, 1 , , , 9 3 , , 2 9 , 9l . . 

-<j)(j) rr + -((pr) +~(p(pr + -a 2 (j) 2 ]; (13) 



r d 2 o d /-, 2m w ^ 2 2 ^ n^t ?ai?e 

{tt - a ~ —T + w^e l + 2 ^o (1 - ' )]}("T~9 - « - -j-)G = 

dr z rdr dr z rdr r 



3 h 3 5- + Q (p(pQ T - (prrPOr H • (14) 

and the boundary conditions 

D r (l) + \[<t>rr{l) + KAV] ~ 3^0 = 0, 

Er (l) + 1^(1) - ^ = 0, £(1) + ^(1) = - J; (15) 

G r (l) + = Kn[4K - ^(1) + MH], G(l) = = . ( 16) 

Here, Si = Si(r) = <fxf>* r - <j) r r<f>* , S2 = ^(f) = §<K ~ 4>r4>* ■ The equations for E* are conjugate 
to those for E. We can use these equations to determine the solution up to 0{e 2 ). The boundary 
conditions are not linearized. Equations (14-16) account for the effect of satisfying the velocity- 
slip condition at the wavy interface, rather at the mean position of the interface. 
To illustrate the nature of the solution, particularly with respect to the effects of nonlinearities 
in convective acceleration and boundary conditions, we shall consider the important case of 
pumping in the absence of the zeroth-order pressure- gradient; i.e., when (dp/dz)o = 0. However, 
a simple superposition of a pressure-driven flow corresponding to equation (7) with Kq ~ 0(e 2 ) 
is permissible. 
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2.1 Originally Quiescent Environment 

To simplify the approach and obtain preliminary analytical solutions of above complicated 
equations and boundary conditions, we only consider the case in which (dp/dz)o vanishes or 
Kq = ipo = 0. Kn is presumed to be comparable with the order of magnitude 0(e) and will be 
rescaled firstly (like Kn) and then adopted by the same representation in the following. Hence 
equations (9-10) become 

(T^-T-- a2 )(T^-T--« 2 )^ = °> a 2 = a 2 -iaRe, (17) 

dr z dr dr z dr 

&•(!) = 0, m = -\, (18) 

V 2 V 2 ^o = 0, 00,(1) = 0; (19) 

together with the condition that the velocity, i.e., {cf> r /r)\ r= Q remains finite along the axis of the 
tube (at r = 0). After lengthy algebraic manipulations, we obtain 

Or = Q r(l-r 2 ); = Arh{ar) + Brh(ar), (20) 

where Q$ is an integration constant; A = Ao/det, B = Bo/det; Aq = aIo(a)/2, -Bo = —aIo(a)/2; 
and 



del = aIo(a)Ii(a) — aIo(a)Ii(a). 
Meanwhile, equation (12) becomes 



2D rrr 3D rr 3D r dL(r) L(r) _ - _ D rr D r 

1 9 q - — j 1 — <-H r J> M r J = —Urrrn n"- 

ry ry£t fy • I CVf* IT IT 7" 



where S(r) denotes the right-hand-side (RHS) term of equation (12). Because of equations 
(12,20), equation (21) can be directly solved, if we denote S(r) = —a 2 Re 2 S(r), so that we have 

n/ . f r -S(t)t 3 dt r 2 f r -, ^ , , x , r 4 r Sit) , log(r)r 2 f r - , 

b 4 + b 3 r 4 + b 2 r 2 \og(r) + b ir 2 , (22) 

where bi,i = 1,2,3,4 are integration constants. To obtain a simple solution which relates to 
the mean flow so long as only terms of 0(e 2 ) are concerned, we see that if every term in the 
z- momentum equation is averaged over an interval of time equal to the period of oscillation [9], 
we obtain for our solution as given by above equations the mean pressure gradient 



Op dp J2,dp. . nu 3,_ C r , 0Orr fori ,2, 1 / n , D rr 



^) + ^[M>*„ - </>*</>„ ~ -{<t><t>* ~ <t>*<t>r)\} + 0(e 3 ). (23) 

ry ry ry 
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Thus, as far as the mean flow is concerned D is the only term which participates in the solution 
as long as only terms of O(e) are retained. We have no need to consider E when considering the 
mean free flow. Actually, with equations (12,21) and the rearrangement of both sides, we have 



dp _ 4Q ,dp _ Pq 

" Tte ' W 2 " W (24) 

where Pq is an integration constant by considering the recombination of equations (12,21) and 
the integration hereafter (cf. equation (23)). 
Now, from equation (14), we have 

D r (l) = -^„.(l) + #r(l)] ( 25 ) 

here, we denote ( = — [</> rr (l) + (f)* r (l)]/2. 

From equation (22) and differentiating D once, we obtain 

D r {r) = (oi + a 2 )r + 2a 2 r In r + a r 3 + G(r), G(r) = -a 2 Re 2 G(r), (26) 

where a 2 must be zero as the axial velocity which is proportional to D r /r remains finite at r = 0. 
Besides, from the expression of D(r), we know that G{r)/r = at r = 0. 

From equations (22-23,25), we know that Pq = — 4ciq. ao or Pq is then proportional to the 
time-averaged pressure-gradient (on the axis) accompanying the peristaltic motion and now 

D r {r) = [C - G(l)]r + a [r 3 - r] + G(r). (27) 

With ao or Po specified, the solution for the mean axial velocity (averaged over time) is 

U EE U(T) = -C^r _ = £[ _Q o(1 _ ,2)] + e 2 { _ [( _ _ G(r) + _ r2)} = 



e[-Q (l " r 2 )] + e 2 {-«o + G(l) - ^ + a (l - r 2 )}. (28) 

r 

In practical applications we must determine ao from considerations of conditions at the ends of 
the channel. As usual, a critical reflux condition can be defined as the condition where there is 
zero velocity at the center of the tube. Using equations (22-23,27), we have 



(§^ = -^[-^ + £(1)], (29) 

where Rq becomes Q (cf. Ref. 19) as Kn=0 and is due to the nonlinear effect of the boundary 
condition. 

3 Results and Discussion 



Our numerical calculations confirm that the mean streamwise velocity distribution (averaged 
over time) due to the induced motion by the wavy elastic interface in the case of free (vacuum) 
pumping is dominated by Rq (or Kn) and the parabolic distribution — Pq(1 — r 2 ). Rq which 



8 



defines the boundary value of D r has its origin in the y-gradient of the first-order streamwise 
velocity distribution. Note that the Reynolds number here is based on the wave speed. The 
physical trend herein is also the same as those reported before 13-14 ' 20 for the slip-flow effects. 
The slip produces decoupling with the inertia of the wavy interface. 

Now, let us define a critical reflux condition as one for which the mean velocity U (r) equals 
to zero at the center-line r = 0. With the equation of U, we have Pq ct = Re(dp/dz)2 = 
— 4[a 2 Re 2 G(l) + Rq], which means the critical reflux condition is reached when Pq has above 
value. Pumping against a positive forcing greater than the critical value would result in a back- 
ward transport (reflux) in the central region of the stream. This critical value depends on a, 
Re, and Kn. There will be no reflux if the pressure gradient is smaller than this Pq. Thus, for 
some Pq values less than Pq ct , the superflow will keep moving forward. On the contrary, parts 
of the flow will move backward if Pq > Pq ct . . 

As reported in Refs. 8 or 17, the rather small critical velocity (< 20/um/s) observed shows 
an apparent dissipation or attenuation of the superflow. Thus, we present some of the val- 
ues of Po(a,Re;Kn = 0,0.15) corresponding to freezed or zero-volume-flow-rate states (cf. 
Chu in Ref. 15) (/„* U(r)dr = 0) in Table 1 where the wave number (a) has the range : 
0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8 and 0.9; the Reynolds number (Re)= 0.1,1,10,25,50 and 100. 
These freezed states might be similar to that reported in Ref. 16 (the nonclassical rotational 
inertia signal (NCRI) is not a universal property of solid 4 He but can be eliminated through an 
annealing of the solid helium sample) . 

We observe that as Kn increases from zero to 0.15, the critical Pq decreases significantly (cf. 
Fig. 2). The quantum slip parameter 13 which is represented in terms of Kn (via the mean free 
path 13 ' 14 ) is temperature dependent in essence but is only relaxed in the boundaries (weakly 
compressible cases here). This kind of boundary-dominated flows resemble those reported in 
Ref. 18 by Sasaki et al. For the same Kn, once Re is larger than 10, critical reflux values Po 
drop rapidly and the wave-modulation effect (due to a) appears. The latter observation might 
be interpreted as the strong coupling between the interface and the inertia of the streaming 
superflow. The illustration of the velocity fields for those zero-flux (zero-volume-flow-rate) or 
freezed states are shown in Figure 3. There is one wave number : a = 0.8 and the Reynolds 
number is 10. Both no-slip and slip (Kn=0.15) cases are presented. 

Some remarks could be made about these states : the transport being freezed in the time- 
averaged sense for specific dissipations (in terms of Reynolds number which is the ratio of 
wave-inertia and viscous shearing effects) and wave numbers (due to the wavy interface or other 
fluctuations) for either no-slip and slip cases. These resemble those reported quite recently 16-18 . 
Meanwhile, the time-averaged transport induced by the wavy interface is proportional to the 
square of the amplitude ratio (although the small amplitude waves being presumed), as can be 
seen in Eqn. (11), which is qualitatively the same as that presented in Ref. 11 for analogous 
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interfacial problems. 

In brief summary, the entrained transport (either postive or negative and there is possibility : 
freezing) due to the wavy interface is mainly tuned by the Pq for fixed Re. Meanwhile, Pq C7 , 
depends strongly on the Knudsen number (Kn, a rarefaction measure) instead of Re or a. These 
results (cf. Table 1 and Fig. 2) might explain why there are rather small critical velocities for 
superflows in the temperature range (similar to the effect represented by the quantum slip pa- 
rameter : Kn here) where a supersolidity is observed 8,17 or the disappearance of supersolidity 
through an annealing of the solid helium sample 16 . We shall investigate much more complicated 
problems 21 ' 22 in the future. 
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Fig. 2 Demonstration of Kn, Re and a effects on the Ho (zero-flux states). 

Re is the Reynolds number (the ratio of the wave-inertia and viscous shearing dissipation). 

a is the wave number and Kn is the Knudsen number (a rarefaction measure). 
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r : distance from tube axis 



Fig. 3 Demonstration of the zero- flux states : the mean velocity field U(r) for 
wave numbers a = 0.8. The Reynolds number is 10. Kn is the rarefaction measure 
(the mean free path of the particles divided by the characteristic length) . 
The integration of U (r) w.r.t. r for these velocity fields gives zero volume flow rate. 



Table 1: Zero-flux or freezed states values (Pq) for a cylindrical wavy interface. 







Re 












Kn 


OL 


0.1 


1 


10 


25 


50 


100 





0.1 


24.2679 


24.2720 


24.2545 


24.1846 


23 9435 

_ '.7 • * J A '.J 'J 


23 0755 


0.2 


24.3293 


24.3271 


24.2753 


24.0009 


23.1347 


20.8498 


0.4 


24.5712 


24.5731 


24.3596 


23.3741 


21.0765 


17.8457 


0.5 


24.7615 


24.7573 


24.4307 


23.0026 


20.1985 


17.1345 


0.6 


24.9884 


24.9855 


24.5116 


22.6346 


19.5060 


16.7273 


0.8 


25.5830 


25.5745 


24.7686 


22.0000 


18.6721 


16.4353 


0.9 


25.9523 


25.9414 


24.9380 


21.8046 


18.4807 


16.4546 


0.15 


0.1 


11.9953 


11.9945 


11.9824 


11.9298 


11.7435 


11.1254 


0.2 


12.0359 


12.0359 


11.9906 


11.7865 


11.1632 


9.7550 


0.4 


12.2082 


12.2066 


12.0421 


11.3259 


9.8944 


8.2432 


0.5 


12.3385 


12.3366 


12.0825 


11.0743 


9.4247 


7.9058 


0.6 


12.4980 


12.4952 


12.1336 


10.8404 


9.0827 


7.6943 


0.8 


12.9223 


12.9152 


12.2967 


10.4917 


8.7025 


7.4861 


0.9 


13.1854 


13.1763 


12.4088 


10.3877 


8.6208 


7.4413 



